Spectroscopy of the fractional vortex eigenfrequency in a long Josephson 0-k junction 
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Fractional Josephson vortices carry a magnetic flux <J>, which is a fraction of the magnetic flux quantum 
•to ~ 2.07 x 10~ 15 Wb. Their properties are very different from the properties of the usual integer fluxons. In 
particular, fractional vortices are pinned and have an oscillation eigenfrequency which is expected to be within 
the Josephson plasma gap. Using microwave spectroscopy, we investigate the dependence of the eigenfrequency 
of a fractional Josephson vortex on its magnetic flux <I> and on the bias current. The experimental results are in 
good agreement with the theoretical predictions. 
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In contrast to integer Josephson vortices known for more 
than 40 years, fractional Josephson vortices are relatively new 
objects. Fractional vortices exist in the so-called 0-K long 
Josephson junctions (LJJs) (J, [2] consisting of 0-parts, having 
the usual current-phase relation (CPR) j s = / c sin(/j), and K- 
parts, having the CPR shifted by K, i.e. j s = j c sin(/j — k). 
Here j s and j c are the supercurrent and the critical cur- 
rent densities and p is the Josephson phase. Often K is 
equal to ?t, and one speaks about 0-7C LJJs, which can be 
fabricated using various LJJ technologies, e.g. based on 
d-wave superconductorsOl 0. S 01 or on a ferromagnetic 
barrierOlliHGii. Using an artificial trick with two tiny 
current injectors, one can also study 0-K LJJs with arbitrary 
K (and fractional vortices), which, in addition, can be tuned 
during experimentffl HHl . 

While an integer vortex (fiuxon) carries one magnetic flux 
quantum <f>o, the fractional vortex carries only a fraction of 
the magnetic flux quantum, i.e. a K-vortex carries the flux 
4> = 4> k/(27c). 

In 0-K LJJs (k > 0) two types of fractional vortices may ex- 
ist: a direct — K vortex, and a complementary K — 2tc vortex| 1 ] . 
The "smaller" vortex corresponds to the ground state of the 
system, while the "bigger" one to the excited state. In the case 
K = 71, both vortices correspond to mirror symmetric semi- 
fluxons with a doubly degenerate ground state. Note, that the 
fiuxon (k = 2k) is an excited state of the system, with a con- 
stant phase being the ground state. Only using topological 
protection, e.g., an annular LJJ, one can trap the fiuxon reli- 
ably. 

In contrast to a fiuxon, the fractional vortex is pinned at the 
0-K boundary. While, being a soliton, the fiuxon may freely 
move along LJJ under the action of various forces (driving, 
friction, magnetic field gradient), the fractional vortex can 
only bend/deform, but it always stays in the vicinity of the 
0-K boundary. When the forces are released, the fractional 
vortex recovers its equilibrium shape. If the damping is small, 
the recovery process is accompanied by decaying oscillations 
around the equilibrium position with the eigenfrequency, as 



predicted in Ref. 

It is crucial to know the eigenfrequencies of a fractional 
vortex. First, any classical device, which uses fractional vor- 
tices, should not operate at frequencies in the vicinity of the 
eigenfrequency (parasitic resonance). Second, an eigenfre- 
quency gives hints about the stability of certain vortex con- 
figurations, namely, a low eigenfrequency is a clear sign that 
the system is close to an instability region, i.e. it may be sen- 
sitive to thermal noise, etc. Third, in the quantum domain 
the eigenfrequency coo defines the attempt frequency of the 
macroscopic quantum tunneling, while, ftcoo defines the en- 
ergy gap between the ground state and the first excited (plas- 
mon) state in the system. 

In this letter we report on the first experimental investi- 
gation of the eigenfrequency of a fractional Josephson vor- 
tex. Using microwave spectroscopy we measure the eigenfre- 
quency as a function of K (or the flux 4> carried by the vor- 
tex) and applied bias current (which deforms the vortex and 
changes its eigenfrequency). 

The eigenfrequency a>o of a single fractional vortex in an 
infinite 0-K LJJ at zero normalized bias current y = I/I c o = 
and vanishing damping a = is 1 12] 
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Here I c o — j c wL is the "intrinsic" critical current, w and L 
are the width and length (circumference) of the LJJ and j c 
is the critical current density. In Eq. Q the prefactor co p o = 
yQojcJ {2kC) is the "intrinsic" zero bias plasma frequency 
related to specific capacitance C of the Josephson barrier and 
j c . For y ^ the analytical expression for a>o(K,y) is un- 
known, but one can approximate it as 
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FIG. 1: Dependences C0o(ft,y) and C0o(37t/2,y) obtained by numeri- 
cally solving the sine-Gordon equation for a 0-K LJJ (symbols) and 
given by Eq. (lines). 



FIG. 3: The dependence 7 e (/i n j) measured at T ~ 4.2K (symbols) 
and corresponding theoretical curve (continuous line). 



where 
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is the normalized critical current of the junction (depinning 
current of the fractional vortex) at given K lll3l Il4l 1 1511 . Ap- 
proximation (|2j differs from the exact numerical solution by 
only few percent [except for y — > Y c (k)], see Fig. \l\ and fol- 
lows the same functional dependence as the plasma frequency 
of a small Josephson i unction Jl6lfl7ll 



«p(y) = « P ov 1 -T ' 
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Thus, the approximation (0 is exact for K = 0. These de- 
pendences serve as a guide for planning and performing the 
experiment. 

For the experiments we used tunnel Nb-A10 x -Nb annular 
LJJ (ALJJ) equipped with two pairs of current injectors as 
shown in Fig. |2] The tunnel ALJJ has low damping and al- 
lows us to perform resonant excitation of the fractional vor- 
tices and perform spectroscopy. The annular geometry, where 
the total topological charge is fixed, prevents flipping of a di- 
rect vortex to a complementary one with emission of a fluxon. 
Even if this happens, upon reset to the zero voltage state, the 
fluxon will be reabsorbed, turning the vortex into the initial 




FIG. 2: Optical image of the investigated sample: ALJJ with two 
pairs of injectors. 



state. Using the injectors we can change the value of K during 
experiment and investigate COo(k). 

The ALJJ has a mean radius R = 30/jm, the width of the in- 
jectors and the distance between them are 5/jm. The Joseph- 
son penetration depth X.t m 43 ^m was estimated taking into 
account the idle repion llal . The normalized length of the 
ALJJ, thus, is txt 4.35. 

To calibrate injectors we have measured the critical current 
I c as a function of 7j n j « K for the left injector pair in Fig. |2] 
This 7 c (/inj) is presented in Fig.[3]and looks like a Fraunhofer 
pattern in accordance with the theory 1 13l ll4lll5ll . The right 
injector pair shows equally good / c (7inj), but was not used in 
this experiment. Knowing 7; n j, the value of K can be calcu- 
lated as K = 27t7 in j/7W n , where 7™j in « 6.92mA is the injector 
current corresponding to the first minimum of 7 c (7; n j)(n|,[3. 

The measurements of the eigenfrequency were performed 
using resonant excitation of the fractional vortex by mi- 
crowaves as follows. We apply a microwave radiation of fixed 
frequency co ex and power P to our 0-K ALJJ with a K-vortex. 
Then we ramp the bias current 7 from zero up to above 7 C and 
observe at which current I\ our ALJJ switches to non-zero 
voltage. Since, according to Eq. (|2ji, the eigenfrequency de- 
creases with the bias current, by ramping 7 one may reach the 
resonant condition coo(K,y) = co ex at some yi = I\/I t ®- The 
fractional vortex will be resonantly excited and will switch 
the ALJJ to the voltage state. In fact, we have been ramping 
the bias current many times, and have measured the escape 
histogram, which represents the probability of switching to 
the voltage state as a function of bias current y, see Fig.|4ji,b. 
These histograms typically have two maxima. The first max- 
imum (at smaller 7) corresponds to the resonant excitation of 
the phase fi(x) (in our case fractional vortex or phase-particle 
from the tilted washboard potential) by microwaves. The sec- 
ond maxima (at larger 7) corresponds to the thermal escape of 
the Josephson phase at 7 approaching 7 c (k). 

Since the accuracy of our measurements is defined by the 
width of the peaks in the escape histograms, special measures 
were taken to suppress electronic noise down to the level that 
the width is due to the bath temperature. We have used a cop- 
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per box directly around our junction to shield electromagnetic 
radiation and a cryoperm shield against static magnetic fields. 
In addition, the measurements were made in a cryostat sur- 
rounded by a three layer /j-metal shield. The whole setup was 
placed in an electromagnetically shielded room. Bias lines 
contained three-stage low pass filters with a 3 dB level cut-off 
frequency at 30kHz (by design): two cold filters are in the 
vicinity of the sample and a warm one at the top of the dip- 
stick. The measurement technique and setup are similar to 
those described in Ref. 20. The bias current was supplied by a 
ramp generator, which starts the ramp at a small negative value 
of current and sends out a digital pulse when the ramp crosses 
zero level. When the ALJJ goes to the non-zero voltage state, 
the voltage detector (embedded in the ramp generator) sends 
out another pulse. By measuring the time between these two 
pulses and knowing the ramp rate, we calculate the switching 
current with high accuracy. The microwave current is induced 
in the bias leads of the ALJJ by means of an antenna placed 
just above the chip inside a copper box and connected to the 
rf generator. A lOdB attenuator is placed right in front of 
the copper box to protect the sample from external noise car- 
ried through the semirigid microwave cable. Furthermore, the 
electronics inside the shielded room is battery powered and 
optically decoupled from the outside[20]. With this setup we 
can obtain a histogram of 10 4 switching events in about 3min. 
The current resolution of this setup is about 30 nA which is 
limited by the noise (jitter) of the current source. 

Before measuring the eigenfrequency of a vortex and com- 
paring it with the theory, we have to measure the plasma fre- 
quency (Opo with reasonable accuracy. To do this we excite 
our system by microwaves and make escape measurements 
at /; n j = 0. In this case the ground state of the system cor- 
responds to a uniform phase fj(x) = arcsin(y) and the lowest 
eigenfrequency is the plasma frequency C0„(y), Eq. (|4}, corre- 
sponding to spatially uniform oscillations (like in a point-like 
junction). For a given microwave frequency co ex the escape 
histogram has two peaks shown in Fig.|^. The first peak at the 
intersection of 03 p (y) and co ex (Fig.0;) corresponds to a reso- 
nant excitation of the phase from the potential well in a tilted 
washboard potentiall 17]. The second peak corresponds to the 
thermal escape of the phase from the potential well when the 
well becomes very 

shallowEEl. 

Its form depends on the 
bias current ramp rate, on the bath temperature T and on the 
electronic noise [2lJ|. Higher T leads to a broader peak shifted 
towards lower currents, which obscures the plasma frequency 
peak especially for low co ex (high y). The electronic noise 
has, to the first order, the same effect as the increase of T. 
Therefore, our setup was optimized as described above to keep 
the broadening caused by the electronic noise below the one 
caused by T. 

The power P of the applied microwaves was kept as low as 
possible, but such that the first peak in the histogram is still 
visible. By increasing P the height of the first peak increases, 
but its position shifts due to the non-linearity of the resonance. 
Accordingly, the height of the second peak decreases and it 
also shifts to lower y. Measurements turned out to be easier 
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FIG. 4: (Color online) Principle of the measurements. The es- 
cape histogram measured at co ex /2lt = 29GHz for (a) K = and (b) 
K = 0.8271 shows two peaks, (c) shows the numerically simulated 
dependences coo(0.82ir, y) and (O p (y) J4}. Applied frequency co ex 
is shown by horizontal line. 



at lower power and large y when the particle is in the shallow 
well and both peaks are close. On the other hand, to resonantly 
excite the particle from the deep well (small y), high power, 
which leads to non-linear effects, is needed. Therefore the 
measurement results for high co ex (low y) are not so accurate 
as for low co ex (high y). 

Escape histograms were measured for different co ex and the 
position of the first peak /i((0 ex ) was, after inverting from 
M<»ex) = M m p) t0 ro p(A) = G>p(y), fitted using the theoret- 
ical dependence @. The fitting gives O3 p o/2n = 42.73GHz 
and a noise free I c q = 961/jA. Note that I c q represents the 
value of j c most accurately. The critical current I c = 7™ ax = 
936/jA measured from the thermal escape peak without mi- 
crowaves represents the maximum supercurrent and has a 
somewhat lower value because of the slightly non-uniform 
current distribution in our ALJJ geometry. The critical cur- 
rent 7 ( IVC = 895 /j A measured from the I—V characteristic or 
from IcQinj) at /mj = (Fig. [3) is even lower. 

An independent numerical study of a microwave driven 
point-like Josephson junction predicts that the position of the 
first peak is somewhat below the plasma resonance 1 22] calcu- 
lated using Eq. For LJJ the situation may be even more 
complicated as the effective ID potential may depend on var- 
ious parameters, e.g., externally applied magnetic field|23] or 
trapped flux. For ALJJ there are no theoretical calculations 
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surements is limited by non-linear effects. For large y (low 
C0 ex ) the measurements are limited by the thermal width of the 
peaks in the escape histogram. To study a>o(K,y) for y — ► 1 
we are going to perform the same measurements at T < 4.2 K. 
Another interesting topic is the spectroscopy of a molecule 
made of two coupled K-vortices, in which eigenfrequency 
splitting is expected. 

We are grateful to the group of Prof. A. Ustinov (espe- 
cially to A. Kemp) from University of Erlangen for sharing 
their experience and numerous suggestions for improving the 
escape measurements setup. This work is supported by the 
Deutsche Forschungsgemeinschaft (projects GO- 1106/1 and 
SFB/TR-21) and by the Eliteforderprogramm of the Landess- 
tiftung Baden- Wiirttemberg. 



FIG. 5: Comparison between the theoretical curves with measure- 
ments of the eigenfrequency. 



or numerical simulations so far. In this work we assume that 
(x) p (y) is well described by Eq. @ even for ALJJ. 

To measure the eigenfrequency of a K-vortex we repeated 
escape measurements for K °< 7; n j > 0. Again, the histogram 
contains two peaks, as shown in Fig. |4j). The first peak 
corresponds to a resonant escape of a fractional vortex at 
COo(k,y) = C0 ex (Fig. 0];). The second peak corresponds to 
the thermally stimulated switching of the junction at the pre- 
critical current. Note, that the critical current of our ALJJ 
depends on K as described by Eq. Q and shown in Fig. [3] We 
have measured many escape histograms for a range of co ex and 
7; n j (i.e. k), and were able to plot coo(k,7). Knowing Q) p o and 
Id) from the plasma frequency measurements we have plotted 
the dependence coo(K,y) in normalized units in Fig. [5] The 
theoretical curves <Mo(y) at different K are shown for compari- 
son. To simulate them, we, first, numerically solved the sine- 
Gordon equation for an ALJJ with injectors of finite width to 
find the static state of the system. Then, to find the eigenfre- 
quency, we have solved the associated eigenvalue problem nu- 
merically and, among all eigenfrequencies, selected the low- 
est positive one (imaginary part of eigenvalue). 

Fig-H]shows that the experimental results are in good agree- 
ment with the theoretical predictions. Comparing the simu- 
lation results for finite injector sizes with the ones for ideal 
(point-like) injectors, we saw that the eigenfrequency in Fig. [5] 
shifts towards larger y (or coo increases for fixed y) as injectors 
get larger. Although we have included the geometrical dimen- 
sions of our injectors in simulations presented in Fig. [5] one 
can see that the experimental points are still somewhat shifted 
towards larger y. This indicates that the effective injector size 
is larger than its geometric size. 

In summary, we have performed spectroscopy of the frac- 
tional vortex eigenfrequency a>o(K,y) in an annular long 
Josephson tunnel junction. The results are in agreement with 
the model. For low values of y the accuracy of our mea- 
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